A percolation transition in Yang-Mills matter at finite number of colours 
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We examine baryonic matter at quark chemical potential of the order of the confinement scale, 
/iq ~ Aqcd- In this regime, quarks are supposed to be confined but baryons are close to the 
"tightly packed limit" where they nearly overlap in configuration space. We show that this system 
wiU exhibit a percolation phase transition when varied in the number of colours N^: at high Nc, 
large distance correlations at quark level are possible even if the quarks are essentially confined. At 
low Nc, this does not happen. We discuss the relevance of this for dense nuclear matter, and argue 
that our results suggest a new "phase transition", varying Nc at constant fiq. 



Long ago, it was suggested to describe the deconfine- 
ment phase transition via percolation [ll-0|- The idea is 
that, at increasing energies, the high parton density will 
make partons of different hadrons overlap. It is natu- 
ral to associate this transition to deconflnement, where a 
quark can propagate throughout the hot medium rather 
than being confined to the hadron size ^qcd 1 f m 
in natural units [8] . Currently, we believe that the transi- 
tion to quark-gluon plasma at low density is a cross-over, 
casting doubt on the relevance of percolation (a phase 
transition, generally of second order) in those conditions. 

In this work we use the percolation picture to study 
a different but related region in the phase diagram, the 
one at low to moderate temperature < T < Aqcd and 
high quark chemical potential fiq Aqcd- Strongly in- 
teracting matter in this regime has recently received a 
considerable amount of interest. Such matter can hope- 
fully be produced in heavy- ion collisions [ol-fl^. and is 
thought to exhibit a rich phenomenology, such as a crit- 
ical point, spinodal instabilities (14j]. separation between 



chiral symmetry and confinement (l5l - [l3. chirally inho- 
mogeneous phases [T9l - [2]| . new phases [22| ■ etc. 

These conjectures are, however, extraordinarily diffi- 
cult to explore quantitatively in a rigorous manner. The 
quark chemical potential /ig is nowhere near the asymp- 
totic freedom limit where perturbative QCD can be used 
Q. It is, however, way too high for existing lattice- 
based approaches, dependent on Hq/T <^ 1, to work 
[23l - [25| . Hence, a simple geometric picture such as per- 
colation might help. Its physical relevance is demon- 
strated in Fig.[Tl schematically showing the regime where 
A*<3 ~ Aqcd- At this chemical potential the density 
is naively expected to be ~ 0(1) Aq^^, that is, one 
baryon per baryonic size. In configuration space, this 
means baryons touch each other, i.e. their quarks are sep- 
arated by a scale not much larger than the confinement 
scale. We therefore expect some weakly-coupled features 
of QCD be present due to the small inter-quark distance. 
Non-perturbative features, on the other hand, should also 
be present since baryons are still confined. The interplay 
of all these features could be very physically interesting. 
Such a setup was recently investigated in [22] using the 
only relevant quantity that can be called "a small param- 
eter": 1/Nc, with Nc number of colours f2^. The idea 



FIG. 1. (colour online) The structure, in configuration space, 
of dense baryonic matter. Comparison of the top and the bot- 
tom panels shows that the percolation picture is applicable: 
when the exchange probability goes above the threshold pc 
the size of the typical cluster diverges. The figure also shows 
a definition of the inter-quark distance y of Eq. [31 



is to keep the quantity A = g'^Nc (where g is the SU {N) 
Yang-Mills coupling constant ) fixed, sending the num- 
ber of colours Nc to infinity, and then expanding in l/Nc- 
It is easy to see that the running of A is qualitatively Nc- 
independent, hence Aqcd ~ Af^- Obviously, such a set- 
up can give at best a qualitative agreement at Nc — 3, 
but it might be enough to understand the phase diagram 
structure of the system. In this setup each quark-quark 
interaction is weak X/Nc) but, due to combinatorics, 
baryons remain strongly-coupled "semiclassical" objects 
[27t . Looking at Fig. [1] with a large- A^c perspective, 
one can immediately see that quarks may be arbitrar- 
ily close together in configuration space (inter-quark dis- 

tance ~ Nc ) , so interactions of quarks between neigh- 
bouring baryons could be weakly coupled. Confinement, 
however, persists until quark-hole screening ~ gluon an- 
tiscreening '^], fJ.q ^ Nc^^Aqcd- 

How can a confined system be, at the same time, 
asymptotically free? The authors of [l^ conjectured that 
excitations below the Fermi surface behave as asymptot- 
ically free quarks, but excitations close to the Fermi sur- 
face are confined. Thus, the entropy of such system is 
^ Nc (unlike N^ for confined matter) , but dynamical ex- 
citations are baryons and mesons and the Polyakov loop 
expectation value Q is ^ 0. This new "quarkyonic" state 
of matter (2^ might be realised in our Nc = 3 world. 

The problem is that it is difficult to model this phase 
theoretically in a rigorous way. While quarkyonic matter 
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has been claimed to be found in the pN JL model [Tl] , one 
can not reliably investigate whether this phase has pres- 
sure ^ Nc in that model, since the pressure in pNJL is 
always somewhat ^ iVp , since the "Polyakov loop field" 
as implemented in [l^ [3 [3 is different from "true" 
confinement. Similarly, it is difficult to see how holo- 
graphic methods [28^1 can verify such a conjecture, since 
there AA^^^^ plays the role of the string coupling con- 
stant gs (A is the compactification scale of the remaining 
5 dimensions in units of Planck length); in a semiclas- 
sical gravity regime, such as_[29j, pressure in a confined 
phase is always ~ [sO, l3l|. Taking these difficulties 
into account, it is not clear how applicable is the large 
Nc limit to our Nc = 3 world. Nuclear matter seems to 
look very different in that world compared to ours, as it 
is a tightly bound crystal of baryons [13, IH, (unless 
what we call large- A^c matter is actually the quarkyonic 
phase). Mean field analysis show a considerable Nc vari- 
ation between TV^ = 3 and 10 > TV^ > 3 [Hi. In [sef 
a physical interpretation of these structures was conjec- 
tured: due to the Pauli exclusion principle, an interplay 
exists between the number of colours Nc and the number 
of neighbours N^. Hence, a "phase transition" exists in 
Nc space when the baryon density is ~ ^qcd i separating 
our world {Nc ^ A^at) from the truly large- A^c world. In 
this work, we aim to apply standard percolation knowl- 
edge to investigate this further. 

The key insight suggesting that interesting structures 
might be lurking in Nc is that 3D bond-percolation ex- 
hibits a phase transition at comparatively low critical link 
probability: for instance, Pc ^ 0.25, 0.18, 0.12 for simple- 
cubic (schematically shown in Fig. [1]), body-centered- 
cubic and hexagonal-close-packed lattice, respectively 
(33 |. Such values suggest that long-distance correlations 
on the quark level could occur even with a somewhat low 
percentage of quarks hopping between baryons, i.e. firmly 
in the confined phase. While below pc the characteristic 
correlation distance ^ (^ cluster size) is ^qcd ' above 
the threshold this quantity explodes to the total system 
size in a comparable amount of time. We leave the mean- 
ing of "correlation" vague, as it can be either a quark hop 
or a gluon exchange; in our context, it implies exchange 
of colour degrees of freedom within a confined tightly 
packed medium. We encode the likelihood of exchange 
between neighbouring baryons in a link probability p, to 
be compared with the percolation threshold Pc in order 
to assess the formation of large-scale structures. 

Two baryons will be correlated if at least two quarks 
are correlated. One has to sum over all possible 
multi-quark configurations, resulting in a strong Nc- 
dependence of p. We determine the latter by calculat- 
ing the probability q — 1 — p that no exchanges happen 
between neighbouring baryons. We define P(i).ij as the 
probability that quarks i and j, respectively in baryons 
A and B, "correlate" (either by fiip or gluon exchange). 
Assuming the quarks inside the nucleon are uncorrelated 
(Fermi motion dominates) , this probability factorises into 
a geometric distribution /yi,_B(x) for quarks to be at a 



certain (vector) position x (see Fig. [T]), and a "squared 
propagator" transition amplitude F[d) for them: 

p=l-(g(i),,,y^=^°; (1) 

= y" /A(xi)dxi j /s(xj)dxj (1 - F(|xi - XjD) . 

We assume a "hard-sphere" distribution for (since 
we keep ^.q fixed, the distance between centres of neigh- 
bouring baryons is always 2Agp^, Fig. [IJ: 

/ab(x) « e(l-AQCD|x-xXT|) (2) 

and a probability of exchange i j based on a range of 
"reasonable" propagators, compatible with confinement 
(fast fall-out in configuration space at distances greater 
than ^ 1 in units of Ag^^) and with the large Nc limit 
of QCD, the interaction is ~ 5^ ^ X/Nc [H, The 
propagators we use are the simple 0-function in config- 
uration space and the momentum-space 0-function used 
in 19], all normalised so their area is Xrx/Nc- In config- 
uration space the transition amplitudes are, respectively, 

A f e(i-^^2£ii) 

ny)-J^{ 2^^^^.fyZ^\ (3) 

Other transition amplitudes, such as a Gaussian distri- 
bution in configuration space, were also tried with no 
significant modifications of the results presented below. 

We note that, due to the fact that 3D percolation 
has a second-order phase transition at a certain pc < \, 
the results we obtain below have some degree of uni- 
versality: as long as the qualitative features of confine- 
ment are observed (the transition amplitude F{y) drops 
sharply above the scale r^, and the hadron density pro- 
file fA,B{'x) has a central plateau of radius ~ ^qcd ^'^d 
a sharply decreasing tail outside), the results we show 
vary quantitatively but not qualitatively, in particular in 
regard to phase transition behavior if p ~ pc for some 
critical Nc- Mathematically, this model is similar to the 
Glauber model, familiar in heavy- ion collisions [1|, with 
the number of colours playing the role of the number of 
participants Npart- Just as in the Glauber model the 
dependence on Npart is universal with respect to cross- 
sections, in this case the dependence on Nc is universal 
with respect to the shape of the transition amplitude. 

The crucial parameter left is a in Eq.[TJ Here, we shall 
consider two limits: one can assume that the link of quark 
i from baryon A to quark j of baryon B does not prejudice 
in any way the possibility of also linking i-k, with quark 
k again in baryon B. This scenario, natural if the link is 
actually realised by a gluon exchange rather than a quark 
fiip, means that a = 2. If, on the other hand, the link 
is given by a quark exchange, then, on a short enough 
time-step, the probability of a quark moving more than 
once is negligible. In this case, a = 1. 

Physically one expects that at A'c ^ 1 gluon exchange 
dominates over quark flip, by combinatorics alone. In- 
deed, one can easily see that a = 1 is in contradiction 
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FIG. 2. (colour online) Contour plot of the critical Nc. for the percolation transition in a hexagonally packed lattice as a function 
of the coupling A and range tt (in Kqcd units). The left panel cissumes a O-function correlation probability in configuration 
space, the right panel assumes a correlation probability based on the propagator used in [ip]. Diagram covers 2 < A^c < 80 



with the Skyrme crystal picture at large Nc- in this pic- 
ture, p{Nc) approaches a constant large-iVc value from 
above: low Nc nuclear matter would be more correlated 
(and hence more strongly bound) than high N^. nuclear 
matter. Comparing strongly coupled Nc — )■ oo nuclear 
matter to the weakly bound nuclear liquid at iVc = 3 
[sBl, this is obviously not right. We therefore assume 
a = 2 henceforward. A closer inspection of the a = 2 
case reveals that it is an approximation of 



p=i - /n/A(x(^')dx(^' //B(x(f))ndx, 
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n(i-nix!^'-xri)) 



(4) 



where we integrate over all quark positions in the two 
baryons and write the probability of no exchanges taking 
place as the product of the individual no-exchanges prob- 
abilities for Nc A-B pairs. Numerical integration shows 
the effect of correlations to be a ^ 3% correction, so that 
the qualitative outcome of the analysis is unaffected. 

In the large Nc limit for the case a = 2, p asymptot- 
ically approaches unity. It is reasonable that this is the 
point where the "dense baryonic matter as a Skyrme crys- 
tal" , theorised in [2^ |23, IH, 111] , is reached. If this is the 
case, however, one should remember that a percolation 
second-order phase transition occurs at apc ^ 1- Hence, 
keeping /ig ^ ^qcd fixed but varying Nc-, the features of 
the Skyrme crystal should manifest not with a continuous 
approach, rather as a second order transition at a not too 
high Nc, whose order parameter can be though to be the 
"giant cluster" density. Below the critical Nc there is lit- 
tle correlation between quarks of different baryons, while 
above this threshold they can correlate, with the distance 
boundary given only by causality. We reiterate that this 

1/2 

is not deconfinement since /i ^ ^qcd ^ Nc ^qcd 
independently of the number of colours, and the frac- 
tion of correlated quarks from different hadrons is still 
~ 0.1 — 0.3 ^ 1 at the percolation transition. Right 
above this transition, therefore, the baryonic wavefunc- 



tion should not be too different from the large- A'c bary- 
onic wavefunction described in [27[ . The correlation dis- 



tance of quarks will however be much larger than the 
baryon size. The features of this new phase are therefore 
similar to those of the quarkyonic matter [23 |. 

Assuming the lowest 3D value Pc — 0.12, appropriate 
for a closely packed hexagonal lattice, the critical number 
of colours is shown in Fig. [5] as a function of A and r^. As 
can be seen, the critical number of colours is significantly 
larger than three for ~ ^q\jd^ ~ 1- However, given 
the roughness of our model, a critical iVc < 3 can not be 
excluded at fiq ~ ^qcd, provided quarks can correlate 
significantly above the confinement length (^ l.SAg^^) 
or the coupling constant is significantly larger than unity 
(one has to remember, however, that these quantities are 
not independent: Aqcd is defined as the scale at which A 
becomes "strong", ~ 1, and it is generally assumed that 
confinement is set by that scale). 

Considering Fig. [5] is a lower limit since pc is at its 
minimum (pc is significantly higher both in a Skyrme 
cubic crystal and in a disordered fiuid), we can say that 
A^c = 3 is disfavoured, although it can not be excluded. 
Changing temperature and /ig should further change the 
critical Nc- Exploring this parameter space, and seeing 
how it relates to the confinement phase transition, is the 
subject of a forthcoming work. 

What are the phenomenological consequences of perco- 
lation? If by "correlation" we mean energy-momentum- 
exchange via quark tunneling between baryons, it is rea- 
sonable that pressure and entropy density ~ Nc above 
the percolation threshold, while below it they stay ~ 
Nc- This is because above the threshold, where inter- 
baryon tunnelling is significant, "typical" excitations of 
the Fermi surface will be superpositions across baryons 
of baryon-localized quark-hole excitations, similar of con- 
duction band electrons in a metal; while the localized 
excitation energy > Aqcd, the superposition makes its 
energy ~ Aqcd even if color degeneracy remains. Thus, 
the degrees of freedom of the system above percolation 
will be delocalized weakly interacting quarks in a lattice 
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of confining potentials, a picture compatible with [20|,|22|. 
Below the threshold, where tunnelling is negligible, exci- 
tations are either color singlets or of energy E ^ Age d , 
suppressed below deconfinenient. 

The only known rigorous way to access these phenom- 
ena quantitatively is the lattice. Current quenched simu- 
lations show that Nc dependence is surprisingly smooth 
[stI . [38| . Our considerations suggest that this will not 
be true at finite chemical potential. While chemical po- 
tentials accessible to current numerical studies are far 
smaller than dense-packing densities |23|j25l| , approaches 
such as the strong coupling expansion [S^ could be used 
to probe the large- iVc dependence at these densities. 

This transition might also be visible with holographic 
techniques beyond the supergravity limit, since finite 
XN~^ corresponds, in gauge/string duality, to gs |28|. 
Thus, percolation will manifest itself as a transition 
from "stringy weak coupling" to "not-so- weak-coupling" . 
While there are hints [40, |41| that percolation is relevant 
for corrections beyond classical gravity, exploring this is 
beyond the scope of this work. The fact that percolation 
appears only as a subleading factor of gs might explain 
why, despite the reasonableness of the argument in [2^ 
for s ^ Nc in the quarkyonic phase, s '-^ in the con- 



fined phase in all semiclassical AdS/CFT setups to date. 

In conclusion, we have used a toy model, with uni- 
versal features, to investigate the close-packed regime 
{fiq ^ ^qcd) of baryons at variable Nc- Our findings 
suggest that, if baryons are kept in this regime but 
is varied, a percolation-type phase transition occurs at 
some critical Nc ^ 0(10), probably but not certainly 
higher than 3. This transition is quite distinct from de- 
confinement as the percentage of quarks propagating on 
superbaryonic distances is quite low. Nevertheless, the 
typical correlation length will be much larger than Ag^^ . 
Further work needs to be done to explore the phenomcno- 
logical consequences of this transition, but our findings 
suggest that applying the Nc ^ oo limit to dense bary- 
onic systems should be done with caution, since a discon- 
tinuity might be present between Nc = 3 and Nc — oo. 
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